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^ I Abstract 

m _ 

<^ \ Comparing Big Bang Nucleosynthesis predictions with the light element abun- 

| dances inferred from observational data, we can obtain the strong constraints on 

(— | ! some neutrino properties, e.g. number of neutrino species, mass, lifetime. Recently 

the deuterium abundances were measured in high red-shift QSO absorption sys- 
tems. It is expected that they are close to the primordial values, however, two 
groups have reported inconsistent values which are different in one order of mag- 
nitude. In this paper we show how we can constrain on r neutrino mass and its 
lifetime in each case when we adopt either high or low deuterium data. We find 
that if 0.01 sec < r Vr < 1 sec and lOMeV < m Vr < 24MeV, the theoretical predic- 
tions agree with the low D/H abundances. On the other hand if we adopt the high 
D/H abundances, we obtain the upper bound of r neutrino mass, m Vr ^z 20MeV. 
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1 Introduction 



It is a very important problem for particle physicists and astrophysicists whether a neu- 
trino has a finite mass or not. Neutrinos are likely massive since it is known that neutrinos 
should have a mass to solve some astrophysical problems such as solar neutrino problem 
(e.g, Among three species of neutrinos, it is presumed that r neutrino is heaviest. 
The present experimental upper bound on the mass of r neutrino is || 



On the other hand, the mass of a stable r neutrino should be less than about lOOeV 
otherwise its cosmic density exceeds the critical density of the universe. Therefore massive 
r neutrino with mass between lOOeV and 24MeV must be unstable, i.e., decays into other 
particles. However the radiative decay of a massive neutrino is stringently constrained 
by laboratory experiments, astrophysics and cosmology and only non-radiative decay is 
allowed. 

Big Bang Nucleosynthesis can strongly constrain on such neutrino properties, in par- 
ticular the number of neutrino species, mass, lifetime, chemical potential and so on, when 
we compare the theoretical prediction of light elements ( 4 He, D, 3 He and 7 Li) with the 
primordial value inferred from the observational data. For example, since neutrino is one 
of the dominant components of the energy density in the universe at the BBN epoch, the 
change of neutrino energy density affects the Hubble expansion rate and alters the freeze 
out time of neutron to proton ratio. Then the predicted light element abundances, espe- 
cially 4 He in this case, are sensitively changed. More precisely we measure the primordial 
component of light element abundances, more strongly we constrain on such neutrino 
properties. 

Recently D abundances were measured in high red-shift QSO absorption systems which 
would not suffer from the effects of the chemical evolution and would have D abundance 
which is closer to the primordial value. However, two groups have reported inconsistent 
D/H values which are different in one order of magnitude. Since we do not have any 
definite models of the galactic and stellar chemical evolutions, it may be premature for 
us to judge which measurements are more correct. 

If the high D abundances, D/H = (1.9±0.5) x 10~ 4 0, in high red-shift QSO absorption 
systems are adopted as the primordial values, since they agree with other light element 
measurements, we can obtain the rigid constraints for some model parameters. On the 
other hand, if the low D abundances, D/H = (3.39 ±0.25) x 1CT 5 j|, [|, are adopted, there 



m VT < 24MeV 
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is a discrepancy between the standard BBN theory and the observational measurements. 
This leads to the similar problem which was pointed out by Hata et al. (1995) || though 
they used the data for D and 3 He in solar neighborhood and a simple chemical evolution 
model to infer the primordial abundances. They found that theoretical predictions agree 
with the observational data if the number of neutrino species N v is less than 2.6. The 
point is that the prediction of the standard (N u = 3) BBN for 4 He is too large if the 
photon-baryon ratio 77 is chosen to fit the D data ( or the prediction of D is too large if rj 
is determined by 4 He data). In this case, we should look for non-standard scenarios. 

There are a couple of ideas proposed to solve the discrepancy. One can reduce the 
abundance of 4 He by assuming non-standard properties of neutrinos (e.g. degenerate 
electron neutrino J7|, massive neutrino ) or can decrease the D abundance by decays of 
exotic particles |§. In particular, it was pointed out || |lOj that a massive and short 
lived r neutrino with mass O(MeV) and lifetime about 0.1 sec. can reduce the effective 
number of neutrino species and hence the abundance of 4 He. In this paper we study the 
effect of a massive and short lived r neutrino on the BBN with a statistical analysis. 
Compared with the recent observational data, especially adopting the low D/H in high z 
QSO absorption systems as a deuterium value, we find that if 0.01 sec ^$ t„ t ^$ 1 sec and 
10 MeV ^5 ^$ 24 MeV the theoretical predictions agree with the all observational data. 
In the case where we choose the high D abundance, we find the upper bound of the r 
neutrino mass, m UT ^20 MeV. 



2 Cosmological Evolution of Massive Neutrino 

We study the BBN effects of the decaying r neutrino by solving a set of Boltzmann 
equations which describe the evolution of momentum distributions of the r neutrino and 
its decay products. A Dirac v T with mass larger than few tens of keV is not allowed since 
it cools SN1987A too fast. Thus, in this paper we only consider a Majorana r neutrino. 
We take into account the three interactions, (1) decay process, (2) pair annihilation and 
creation process, and (3) scattering process. Concerning the decay process, we suppose 
that a massive Majorana r neutrino with mass m Vr and lifetime t Ut , decays into a fi 
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neutrino and a scalar particle:^ 

v T -> ^ M + 0, (2) 

where we assume that </> and z/ M are massless. The scalar particle may be a Majoron [II] or 
a Familon |12j and the decay modes are expected to be detected in the future experiments, 
CLEO III or B factories ]13| . The Boltzmann equation for the distribution /, of i-particle 
(i = u T , u^, 0) is written as ]IJ], |T5|, 

~ir ~ ~T7\Pi^T~ = decay) + d(ann) + CAscatt) (3) 
at a(t) dpi 

where a(t) is the scale factor, p^ is the momentum, and C^s are collision terms which 
correspond to the three processes. 

The collision term C Ut for the decay process v T {jPu T ) — > u^ip,, ) +4>(p</>) is given by || 

C VT {decay) = -r£ +T& (4) 

where and Tj D are the collision terms for the decay and the inverse decay processes, 
respectively. They are given by, 

n (5) 
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The integration is performed between p max = {Jrrfy T + p\ T +p Ur )/2 andp mim = {sjm\ T + p\ T ~ 
p„ T )/2. In the same way, the collision terms C Ull (decay) and C^decay) are given by, 

C^decay) = - rft* (7) 
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2 The t neutrino might decay into a electron neutrino which changes the neutron-proton ratio directly. 
However, in the present paper we do not consider this decay channel since the Boltzmann equation for 
v e is too complicated. It is expected the emissions of v e and i> e increase the neutron-proton ratio and 
tend to enhance the discrepancy between the BBN theory and the observation. 
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where k ^ k' = or 0, the spin factors are g Vr — g u = 2, and — 

For pair annihilation and creation processes v{p v ) + v{p' u ) — > h{Pi) + ^(Pi), 
e, v ei u^), the collision terms are written as 



C(ann) = -— J p' k 2 dp' k (av) k (f k (p k )f k (p' k ) - f eq (p k ) f eq (p' k )) , (10) 

where k = z/ T , and f eq is the equilibrium distribution and (av) k is the differential 
cross sections given by 
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(H„ = ^f E + (12) 



i=e, u e ,u u 



where Gf is Fermi coupling constant, and Ca = 1/2, Cy = 1/2 — 2sin 2 9w (9w- Weinberg 
angle) for and Ca = 1/2, Cy = 1/2 for neutrinos. We do not include the annihilation 
and creation process for the scalar particle <fi assuming that their interaction is very week. 

As for scattering processes, it is not easy to solve the Boltzmann equations including 
the collision terms due to scattering because they are generally much more complicated 
than collision terms due to the annihilation and the decay. Without scattering, the spec- 
trum of the r neutrino is distorted because of the momentum dependence of the annihila- 
tion rate, which leads to overestimation of the density of v T ||]. This effect is significant for 
the neutrino whose mass is larger than about lOMeV since such heavy neutrino becomes 
non-relativistic before decoupling and its density decreases through annihilation processes. 
On the other hand, for neutrinos which are relativistic at decoupling epoch, the momen- 
tum distributions do not change by the annihilation and is the same as the equilibrium one. 
Thus the scattering process is important only for the heavy r neutrino (m UT ^ lOMeV). 
Therefore we consider only the scattering process v T {p) + l(k) — > v T {p') + in the 
collision term of the Boltzmann equation for v T . Using non-relativistic approximation, 
the collision term is simplified as 

C{scatt) = -T f + F b , (13) 
where T f and Tb are forward and backward reaction rates which are given by 

r '^£,^ + ^ m ^w^ hM (14) 
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(15) 
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r(y,y')fi(y'), 



y/x' 2 + y' : 

where x = m UT /T, y = p/T, and y' = p'/T. The function r(y, y') is given by 

/ /\ [a 3, 1. /|2 \ \v-v'\ 

r{y,y) = lA + -\y-y\ + -\y-y\ Je ? 

^ + \{y + y') + \{y + y') 2 ) e _J ^- 



(16) 



At T ^ 0.5 MeV, electrons and positrons begin to annihilate and heat photons without 
heating neutrinos. Then we have to deal with the temperature and cosmic time relation 
carefully. In this situation, we solve a set of Boltzmann equations coupled with the 
Hubble expansion rate and the entropy conservation equation. The Hubble expansion 
rate is obtained by solving the Einstein equation: 
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and the entropy conservation is: 



(17) 
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(18) 



where pi denotes the energy density of the particle "i" (i = e, 7, u e , u T , 0) and P e is the 
electron pressure. 

As initial conditions we take 



fv T — f^ij. ~ feq 

U = o, 



(19) 
(20) 



where we assume that the scalar particle has been never thermalized and its cosmological 
density is negligible. We will discuss the justification of this assumption later. 

3 Time Evolution of N u 

As we have already mentioned, the abundances of 4 He is sensitive to the total energy 
density of the universe during the BBN epoch, T = 10 MeV — 0.1 MeV. The lower the 
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energy density, the more slowly the universe expands, which leads to later decoupling of 
the weak interactions and less neutron-proton ratio. Then the abundance of 4 He decreases 
since almost all neutrons are translated into 4 He. Since the other light elements are 
insensitive to the cosmic expansion rate, it is expected that the BBN predictions may 
agree with observational data if the energy density of the universe becomes lower than 
the standard one by the effect of heavy r neutrino decay. 

When we discuss the energy density of massive decaying r neutrino, it is useful to 
introduce the effective number of neutrino species N v defined by 

N v = P-.+Pr+P-r+P* . (21) 
p v [m v = 0) 

The denominator is the energy density of a massless neutrino and N v is a constant ( = 3 
) in the standard model. 

In Fig. 1 we show the time evolutions of N v for (t Vt ( sec), m^ T (MeV)) = (10~ 3 , 1), (1, 10) 
and (0.1, 14). In order to understand their behaviors, it is convenient to introduce the three 
temperature scales, (1) the epoch when the decaying particle becomes non-relativistic 
T NR , (2) the decaying particle's decoupling epoch T D , and (3) the epoch when the decay 
occurs T decay . We define them by, 

T NR = (22) 
T(T D ) = H(T D ), (23) 
T decay = T(t = r eff ), (24) 

where H{t) is the Hubble expansion rate, and r e ff is the effective lifetime of the r neutrino. 
Notice that the effective lifetime of the relativistic neutrino becomes longer than the proper 
lifetime because of the time dilation ( r e ff ~ T^m 2 VT + V 2 v r l m v T \ 

In the case of t Vt = 10 -3 sec and m VT — 1 MeV, the relation of the three time scales 
is given by T decay > T D > T NR . At temperatures between T decay and T D the decay, inverse 
decay, annihilation and pair creation processes tend to establish the thermal equilibrium 
among three species of neutrinos and the scalar particle and N v becomes close to 3.6. 
Then, after decoupling, N v increases due to the mass effect until T^r at which the inverse 
decay becomes suppressed and v T begins to decay efficiently On the other hand, in the 
case of t Vt = 1 sec and m UT = 10 MeV, T NR > T D > T decay . Between T NR and T D the 
energy density of r neutrino decrease due to the Boltzmann suppression. After decoupling, 
the mass effect increases N u , and at T ~ T decay v T begins to decay. Thus, if Tjj ~ T decay , 
the increases of N v due to the mass effect is not significant and becomes less than 
three as shown in Fig. 1 for t Vt = 10 _1 sec and m VT = 14 MeV. 
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4 BBN with Massive Decaying r Neutrino 



In a previous section, we have discussed the time evolution of the energy density in the 
universe, and have shown that r neutrino with O(10 MeV) mass and short lifetime (~ 0.1 
sec) can reduce the effective number of neutrino species N v . In order to calculate the 
abundances of D, 4 He, 3 He, and 7 Li, we incorporate the change of the effective number 



of neutrino species N„ into the big bang nucleosynthesis code |L0 . 

In the case of massive and unstable tau neutrino, the effective number of neutrino 
species is a time-dependent variable. Since only 4 He abundance is sensitive to N„, if the 
predicted 4 He abundance is expressed by the corresponding (constant) number of neutrino 
species in the BBN calculation without massive neutrinos, it helps us to understand the 
results of our calculation intuitively. Thus, we introduce the equivalent number of neutrino 
species N v here which is defined as 

Y p (N ul m Ur = 0,t„ t = oo) = Y P (N V = 3, m Vr , t Ut ). (25) 

In Fig. 2 we show the contour plots of the equivalent number of neutrino species N v in 
the parameter space, t Ut —m UT . In this plot we expect that the predicted Y p decreases if 
t Ut ^ 10 sec and m Vr > 10 MeV. 

In the numerical calculation we make use of the Monte Carlo procedure in order 
to include the experimental uncertainties of the nuclear reaction rates and the neutron 
lifetime [fT?! [HI - Since the theoretical la errors are almost independent of N v for N v = 
2 — 4, we adopt the la errors estimated by the standard BBN (N v = 3) calculation. 

The calculated abundances are compared with the primordial values inferred form the 
observational data. Concerning the 4 He, we adopt the value in the extra galactic HII 
region and it is given by 



Y p = 0.234 ± 0.002(stat) ± 0.005(syst) (26) 

where Yp is the mass fraction of 4 He and "p" denotes the primordial value. It is believed 
that abundances of 7 Li observed in the population II metal poor halo stars in our galaxy 



represent the primordial values. We adopt the recent observational data p0 



( 7 Li/H) p = (1.73 ± 0.05(stat) ± 0.20(syst)) x 10" 10 . (27) 

For a deuterium abundance, we have some types of the measurements which are the 
local ISM or the solar system data and the abundances observed in high red-shift QSO 
absorption systems. The abundances inferred from the local ISM are (D/tfyisM = (1-6 ± 
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0.2) x 10 ||21||. In the solar system the observed deuterium abundances are (D/H) G 







2.57 ±0.92) x 10 5 p2 |. This pre-solar deuterium abundance is obtained by subtracting 



3 He/H value observed in carbonaceous chondorites from (D + 3 He)/H observed in the 
solar winds. From above abundances, we must extrapolate the true primordial deuterium 
abundance assuming appropriate chemical evolution models. 

On the other hand, The observations in high red-shift QSO absorption systems are 
expected that they are nearly the primordial values because such Lyman a clouds would be 
unprocessed. In the last three years, more papers concerning these Lyman limit system 
observations are published. However, two groups have measured different D/H values 
which disagree in a order of magnitude. The higher values measured by some persons 
are 0, 

(D/H) qso = (1-9 ±0.5) x 10" 4 . (28) 
The lower values measured by Buries and Tytler. || are, 

(D/H) qso = (3.39 ± 0.25) x lO" 5 . (29) 

Some chemical evolution models which predict only modest destruction of deuterium 
during the galaxy evolution have been suggested |2ij|. Adopting such a recent model 



| 24fl which reports that the deuterium destruction is limited to a factor 3 or less, the above 
low deuterium values in high red-shift QSO absorption systems would be consistent with 
the ISM and solar system data. However, there are also the other models in which pri- 
mordial deuterium can be significantly destroyed into a order of magnitude lower without 
over-producing the 3 He |[25|| . The latter models might agree with the high deuterium 
abundance in high z QSO absorption systems. Considering the circumstances mentioned 
above, we so far have no conclusive model of the galactic chemical evolution. On the 
observational side, it is premature to determine which deuterium values in high red-shift 
QSO absorption systems are reliable. Therefore in this paper, we adopt both the high 
and low deuterium data in high z QSO absorption systems and for each deuterium value 
we constrain on the r neutrino mass and its lifetime. 

In order to obtain the confidence level of the model parameters, mass m UT , lifetime 
r„ T , and baryon to photon ratio 77 in the universe, we utilize the Maximum Likelihood 
analysis. The likelihood function of theoretical abundances is taken to be gaussian which 
has the one sigma variance calculated by the Monte Carlo procedure. The likelihood 
function of observational data is taken to be gaussian for the statistical error, and to be 
top hat distribution for the systematic error (see Ref. [0]). Also for the r neutrino mass, 
the experimental bound Eq. ([I]) is incorporated as a gaussian distribution. 
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The obtained two dimensional contours in the parameter space, r„ T —m UT , are shown 
in Fig. 3 for using the low (D/H)q SO , and in Fig. 4 for using the high (D/H) qso . In 
Fig. 3 the contour with same likelihood as the standard BBN model (N^ = 3) is shown 
by the dotted line. Notice that if 0.01 sec < r„ T < 1 sec and 10 MeV < m„ T < 24 MeV, the 
theoretical predictions give a much better fit to the observational data than the standard 
BBN. The best fit values are estimated as 

m Vr = 15 tl 1 MeV (68% C.L.), (30) 

logiofe) = -1-0 ±i:S (68%C.L.), (31) 

for the low D abundances in the QSO absorption systems. The filled square in Fig. 3 
represents the best fit point. From Fig. 4 we can obtain the upper limit for the r neutrino 

mass m„ r , 

m VT < 20 MeV (68% C.L.), (32) 

for high D abundances in the QSO absorption systems. 



5 Conclusion and Discussion 

In this paper we have shown the light element abundances which are predicted by BBN 
theory with r neutrino whose mass is ~ lOMeV and lifetime is ~ 0.1 sec. agree with the 
primordial values inferred by observational data even if we adopt the low D/H in high z 
QSO absorption systems. On the other hand, if we adopt the high D/H, we can obtain 
the strong constraints on r neutrino mass m UT < 20MeV. 

Finally we discuss the annihilation rate for the scalar particle 0. For this purpose we 
assume that is a Majoron-like particle which has the interaction to u T as 

Lint = -d^Ur'ys'fUr, (33) 
V 

where v is the characteristic energy scale. In the relativistic limit the ratio of the annihi- 
lation rate for + — > v T + v T to Hubble expansion rate is given by 

I^^.^^'^^g (T» mJ . (34) 

When the cosmic temperature T is less than m Ur , the ratio Y/H oc T 2 . Thus Y/H takes 
the maximum value at T ~ m VT : 

4 
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The condition that the scalar particle is not in thermal equilibrium is given by 



w>7xl0 4 MeV (36) 

On the other hand, the lifetime of the r neutrino is related to v by 

16irv 2 , , 

(37) 



2 - 



where U T)I1 is the mixing angle. Therefore the condition ( |36"D is written as 

V>4x 10" 10 sec \U T J- 2 . (38) 

From Eq. ( j3~8"D the scalar particles are not in thermal equilibrium if t Ut > lxlO- 5 sec and 
\U T J> 10' 2 . Thus our assumption that the scalar particle has not been thermalized is 
justified when the mixing angle is large. 

We wish to thank T. Asaka for useful discussions for the Majorana neutrino. This 
work has been supported by the Grant-in-Aid for COE Research (07CEQ2002) of the 
Ministry of Education, Science, and Culture in Japan. 
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Figure caption 

Figurel: Evolution of the effective number of neutrino species N v with the cosmic tem- 
perature T for (r i , T (sec),m ! , T (MeV)) = (KT 3 , 1), (1, 10) and (0.1, 14). 



Figure2: Contours of the equivalent number of neutrino species N v in the parameter 
space, t Ut — m VT . N v is defined by Eq. (pop. In this plot we expect that the predicted 
Y p decreases if t Vt < 10 sec and m Vr ^ 10 MeV. By the definition, N v does not depend on 
the baryon to photon ratio r\ at all. 



Figure3: Contours of the confidence levels in the parameter space, t Vt — m Vr . In this 
plot the lower deuterium abundance obtained in high red-shift QSO absorption systems 
is adopted as the primordial value. The dashed line denotes 90% C.L., and the solid line 
denotes 68% C.L.. The contour with same likelihood as standard BBN model (N„ = 3) 
is shown by the dotted line. The filled square represents the best fit point. 



Figure4: Contours of the confidence levels in the parameter space, t Vt — m VT . In this 
plot the higher deuterium abundance obtained in high red-shift QSO absorption systems 
is adopted as the primordial value. The dashed line denotes 90% C.L., and the solid line 
denotes 68% C.L.. The contour with same likelihood as standard BBN model (N„ = 3) 
is shown by the dotted line. 
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